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Abstract 

We present a new measurement of the branching ratio R of the decay Kl — > ir^e^v, 
denoted as K e s, relative to all charged Kl decays with two tracks, based on data 
taken with the NA48 detector at the CERN SPS. We measure R = 0.4978 ± 0.0035. 
From this we derive the K e ^ branching fraction and the weak coupling parameter 
\V US \ in the CKM matrix. We obtain \V us \f+(0) = 0.2146 ± 0.0016, where /+(0) is 
the vector form factor in the decay. 
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1 Introduction 



The unitary condition for the first row of the CKM quark mixing matrix 
is at present fulfilled only at the 10% C. L. [1]. This has renewed interest 
in the measurement of the coupling constant V us for strangeness-changing 
weak transitions. The most precise information on V us comes from the decay 
Kl — > 7r ± e =F z/, which is a vector transition, and therefore is protected from 
SU(3) breaking effects by the Ademollo-Gatto theorem [2]. We present here a 
new measurement with improved experimental precision. 



2 Apparatus 

The experiment was performed using the NA48 detector in a beam of long- 
lived neutral kaons produced at the 450 GeV proton synchrotron SPS at 
CERN. The neutral Kl beam was derived at an angle of 2.4mrad from an 
extracted proton beam hitting a beryllium target. The decay region starts at 
the exit face of the last of three collimators 126 m downstream of the target. 
The experiment was originally designed and used for the precision measure- 
ment of direct CP violation in kaon decays [3]. We report here on a study of 
semileptonic decays, for which data were taken in a pure Kl beam in Septem- 
ber 1999. The main elements of the detector relevant for this exposure are the 
following: 

The magnetic spectrometer consists of four drift chambers (DCH), each with 
8 planes of sense wires oriented along four projections, each one rotated by 
45 degrees with respect to the previous one. The spatial resolution achieved 
per projection is 100 /im, and the time resolution for an event is 0.7 ns. The 
volume between the chambers is filled with helium near atmospheric pressure. 
The spectrometer magnet is a dipole with a field integral of 0.883 Tm and is 
placed after the first two chambers. The distance between the first and the last 
chamber is 21.8 meters. The spectrometer is designed to measure the momenta 
of the charged particles with high precision - the momentum resolution is given 
by 

a(p)/p = 0.48% ©0.009 ■ p% (1) 
where p is in GeV/c. 

The hodoscope is placed downstream from the last drift chamber. It consists 

GZ 616.360/2-IV GZ 616. 363/2- VIII, and by the Austrian Science Foundation under 
contract P08929-PHY. 
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of two planes of scintillators segmented in horizontal and vertical strips and 
arranged in four quadrants. The signals are used for a fast coincidence of two 
charged particles in the trigger. The time resolution from the hodoscope is 
~ 200 ps per track. 

The electromagnetic calorimeter (Lkr) is a quasi-homogeneous calorimeter 
based on liquid krypton, with tower readout. The 13212 readout cells have 
cross sections of ~ 2 x 2 cm 2 . The electrodes extend from the front to the 
back of the detector in a small angle accordion geometry. The Lkr calorime- 
ter measures the e ± and 7 energies by summing the ionization from their 
electromagnetic showers. The energy resolution is: 



where E is in GeV. 

Charged decays were triggered with a two-level trigger system. The trigger 
requirements were two charged particles in the scintillator hodoscope or in 
the drift chambers coming from a vertex in the decay region. 

A more detailed description of the NA48 setup can be found elsewhere [3]. 



3 Data analysis 

3.1 Analysis strategy and events selection 

The basic quantity measured in this experiment is the ratio R of decay rates 
of K e3 decays relative to all decays with two charged particles in the final 
state, mainly neu, njiu (called K^), 7r + 7r~7r° (called K 3n ), 7r + 7r~ (called K 2n ) 
and 37T° with Dalitz decay of one 7r°, denoted as 7r°7r°ee7 or 7r°7r 7r^,. Since 
the neutral decay modes to 3tt°, 2tt° and 77 have been measured, and the 
correction for events with four tracks B(4T) is small, the sum of branching 
ratios of all K L decay modes with two charged tracks B(2T) is experimentally 
known [1] 



a(E)/E = 3.2 %/VE © 9.0 %/E © 0.42 % 



(2) 



B(2T) = 1 



= 1 




(3) 
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Since contributions from K$ meson or A hyperon decays and from rare Kl 
decays are negligible, less than 2 x 10~ 5 , we obtain the K e3 branching ratio 
using B{2T): 



B(e3) 



x B(2T) . 



(4) 



r(K L 



all) T(K L -> all 2-track) 



In this experiment, we therefore measure the ratio of K e3 events to all 2-track 
events N 2 t divided by their acceptances a e or a 2 r respectively: 



Both numbers, N e and N 2 t are extracted from the same sample of about 80 
million recorded 2-track events. These were reconstructed and subjected to 
offline filtering. 

In the basic selection, two tracks were required with opposite charge and 
a distance of closest approach below 3 cm. The vertex fiducial volume was 
defined to be between 8 m and 33 m from the final collimator, and within 3 cm 
of the beam-line. Events with high hit multiplicity were rejected by requiring 
that no overflow condition occurred in the drift chambers. An overflow is 
generated if more than seven hits in a plane were recorded within 100 ns. 
These cuts were passed by 48.795 million events. 

Events were rejected if the time difference between the tracks exceeded 6 ns. 
Both tracks were required to be inside the detector acceptance and within 
the momentum interval lOGeV/c to 120GeV/c. In order to allow a clear 
separation of pion and electron showers, we required the distance between 
the entry points of the two tracks at the front face of the electromagnetic 
calorimeter to be larger than 25 cm. 

The last selection criterion was applied to a measure of the kaon momentum, 
to avoid the region below 50 GeV/c which is simulated inadequately. We used 
the sum of the moduli of the two momenta P = P\ + P 2 for all decays. As a 
result 12.592 million events with P > 60 GeV/c remained. For the denomina- 
tor, no identification of individual decay modes was applied, and the average 
acceptance a 2 r applies to the requirements listed up to this point. 

For the numerator A^, the K e3 signal was selected by a single additional cri- 
terion that at least one track should be consistent with an electron. This was 
done by requiring that the ratio E/p exceed 0.93, where E is the measured 
energy in the calorimeter and p is the measured momentum in the magnetic 
spectrometer. 6.759 million events were accepted. The quantity E/p is shown 
in Fig. 1 for all tracks of these K e3 events. 



R = 



N e /a e 



(5) 
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0.2 0.4 0.6 0.8 1 1.2 

E/p: energy deposited in Lkr / track momentum 

Fig. 1. The ratio of calorimetric energy E over the momentum p 
for the tracks of all selected K e ^ events. 



3.2 Corrections for electron identification 



The number of K e ^ events was corrected for the inefficiency of the electron 
identification (electrons with E/p < 0.93) and background coming from 
and decays (pions with E/p > 0.93). Both inefficiency and background 
were measured from the data. 

For the background determination a sample of events was selected having one 
track with E/p > 1.0, clearly classifying it as an electron. The background 
probability for pions W(7r —> e) was then determined from the E/p spectrum 
of the other (i.e. pion) track (see Fig. 2) to be 

W{jx -»• e) = (0.576 ± 0.005(stat.)) %. 



As a cross check the probability was also derived from the E/p spectrum 
of K 37T events, giving a consistent result within errors. Background from the 
decay Kl — » 7r°7r 7r£, was completely removed by the cut on P. 

The electron ID inefficiency W(e — > tt) was determined in a similar way (see 
Fig. 3) by requiring one track with E/p < 0.7, classifying it as a pion. The 
E/p distribution for the other track then consists mainly of electrons, with a 
small contribution from pions, especially below 0.7. We subtracted this pion 
component by using the previously determined pion distribution, normalized 
in the range 0.2 < E/p < 0.6. From this we then obtained the probability for 
losing an electron by the condition E/p > 0.93: 

W(e -> tt) = (0.487 ± 0.004(stat.)) %. 
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0.2 0.4 0.6 0.8 1 1.2 

Ke3 events: E/p of pion tracks 



Fig. 2. Quantity E/p for pion tracks. The sample was selected by 
the requirement E/p > 1.0 for the other (i.e. electron) track. 




0.2 0.4 0.6 0.8 1 1.2 

Ke3 events: E/p of electron and pion tracks 



Fig. 3. Quantity E/p for electron and for pion tracks. The electron 
spectrum is scaled for better illustration. The dark shaded area 
represents electrons with E/p < 0.93. 

3.3 Monte Carlo Simulation 

To reproduce the detector response, a GEANT-based simulation of the NA48 
apparatus was employed for the five decay modes 7rez/, TTfiu, tt + 7t~7t°, 7t + tt~ 
and 7r°7r 7r^. Radiative corrections were included for the K e3 mode. We used 
the PHOTOS program package [4] to simulate bremsstrahlung, and added 
the calculations from [5] on virtual photons and electrons. Some comparisons 
between data and MC for identified K e3 events are shown in Fig. 4 (z-vertex) 
and Fig. 5 (x- and y-coordinates of the tracks in the first drift chamber). 
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Fig. 4. Longitudinal vertex distribution for K e3 events: data and MC (left) 
and ratio of data over Monte Carlo simulation (right). 

We obtain the individual acceptances shown in Table 1. 



decay mode 


acceptance 


K e3 


0.2599 




0.2849 


K 37r 


0.0975 


K 27r 


0.5229 




0.0001 



Table 1 

Detector acceptances for the charged decay modes. 

The average 2-track acceptance a 2 x was obtained from a weighted mean of 
the individual acceptances a« which depends only on ratios of decay rates 
measured in other experiments: 



d 2 T 



B e a e + B^a^ + B 37r a 3n + B 27r a 2n + B D a D 
B e + B^ + B 37T + B 2w + B D 

e ^ Be a e Be a e Be CL E Be a e I 



B D ' 



(6) 



Here Bi are the branching ratios for the decay channels (i = e : K e3 ; i = fi : K^; 
i = 3iT : 7T + 7r _ 7r ; i = 2tt : 7t + tt~; i = D : ti ^ ^). The acceptance for channel 
i is Oj. For the branching ratios we used a weighted average of the 2004 PDG 
values [1] and the new KTeV measurement [11]. The uncertainty was enlarged 
according to PDG rules for averaging inconsistent data. 
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K e3 : Data/MC 



-80 -60 -40 -20 20 40 60 80 
track Y-coordinate in DCH1 [cm] 



Fig. 5. Transverse positions (horizontal x and vertical y) of tracks in 
the first drift chamber from K e3 events: data and MC (left) and ratio of 
data over Monte Carlo simulation (right). 

BjB e = 0.666 ±0.011 
B 3n /B e = 0.309 ± 0.004 
B 2 n/B e = (4.90 ± 0.14) x 10~ 3 

B D /B e = (1.96 ±0.05) x 10" 2 



(7) 
(8) 
(9) 
(10) 



Varying the constraints given by Eq. (7) to (10) within their errors we get a 
relative variation of the acceptance of 0.16 %, and a^T = 0.2412 ± 0.0004. 



3.4 Systematic uncertainties 



Given the large number of events, the uncertainties of this measurement are 
purely of systematic nature. Simulation shows that most of these systematics 
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will induce a dependence of the result on the lower cut on the sum of the two 
moduli of the two momenta P = Pi + Pi- Since the three decay modes have 
different neutral energy, which is either invisible as a neutrino or not used in 
this analysis, events with a given value of P originate from different average 
kaon energies, so a possible imperfection of the kaon energy spectrum (which 
is fairly well known for energies above 50 GeV) will induce a dependence of 
the result on P. 

In Fig. 6 we show a comparison of the energy spectra for identified K 3n and 
K 2n events, where we can fully reconstruct the energy. Both comparisons show 
a small slope but with opposite signs, demonstrating that the kaon energy 
spectrum in the MC is a good compromise between different decay modes. 

Fig. 7 compares the momenta of the electrons and pions in identified K e3 events 
between data and MC. Fig. 8 shows the same comparison for the sum of track 
momenta in the range between 60 GeV/c and 130 GeV/c, which contains 95 % 





kaon energy [GeV] kaon energy [GeV] 

Fig. 6. Reconstructed kaon energy from and decays, comparison 
between data and MC. Errors are statistical only. 
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of the data. With radiative corrections applied, we still observe a slight slope 
in P of half the size of the slopes of the fully reconstructed events in Fig. 6. 
This is the dominant source of experimental uncertainty, and may be due to 
imperfections of the radiative corrections as well as limitations in the detailed 
event simulation. 



1.2 



1.1 



0.9 



0.8 



0.7 



K e3 : Data/MC 



Electron 



i ... i ... i . i . i 



1.2 



1.1 



20 40 60 80 100 
electron track momentum [GeV/c] 



0.9 



0.8 



0.7 



K e3 : Data/MC 



Pion 



i ... i ... i ... i 



20 40 60 80 100 
pion track momentum [GeV/c] 



Fig. 7. Comparison between data and MC (including radiative correc- 
tions) for the momenta of electrons and pions in identified K e ^ events. 
Errors are statistical only. 



1.1 



1.05 



0.95 



0.9 



K e3 : Data/MC 



Slope: 

(-2.22 ± 0.41) x10" 4 /GeV/c 




lit 



60 70 80 90 100 110 120 130 140 
sum of track momenta [GeV/c] 

Fig. 8. Comparison between data and MC for the sum of track 
momenta. The region between 60GeV/c and 130GeV/c contains 
95 % of the data. Errors are statistical only. 



To get a conservative estimate of this dependence, we varied the lower cut on 
the value of P from 50 GeV/c to 80 GeV/c. This is a large range of variation, 
considering that the analysis used data above 60GeV/c, and that a cut at 
80 GeV/c removes 70 % of the events. The resulting relative uncertainty of the 
ratio R in Eq. (5) is 0.67%. In addition, a second independent analysis was 
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performed using different selection criteria and a different kaon momentum 
spectrum, which was weighted such as to reproduce exactly the kaon momen- 
tum spectrum of K e3 events. The value of R differed from the one in the first 
analysis by 0.2%, well below the estimated systematic uncertainty. 

To estimate the uncertainty coming from the E/p cut to select K e3 events, 
we varied the cut value between E/p > 0.90 and E/p > 0.96. As a result, 
inefficiency and background due to this criteria vary significantly, leading to 
different net corrections of K e3 event numbers (Table 2). Applying these cor- 
rections, however, we get almost the same number of events, thus demon- 
strating the correctness of this selection principle. It appears that with E/p 
> 0.93, both inefficiency and background are very small and nearly cancel. 
The resulting relative uncertainty on R is AR/R = 0.05 %. 





E/p > 0.90 


E/p > 0.93 


E/p > 0.96 


inefficiency [%] 


0.275 


0.487 


1.424 


background [%] 


0.914 


0.576 


0.266 


K e z event number after E/p cut 


6796461 


6759184 


6673114 


net K e s correction 


-42624 


-5705 


77182 


corrected K e3 event number 


6753836 


6753478 


6750296 



Table 2 

Variation of the E/p cut to select K e3 events. 



The data used in this analysis originate from two different triggers (Q2 + 
Ql/20 * 2trk), where Q2 requires two quadrants of the hodoscope counter to 
be hit, while Ql*2trk requires at least one hodoscope quadrant plus two tracks 
from the drift chamber trigger system. Ql is prescaled by a factor of 20. By 
selecting one trigger, the efficiency of the other can be measured, taking into 
account the different downscaling. The trigger efficiencies for 2-track and K e3 
events differ slightly for the Q2 trigger ( (97.38 ± 0.02) % for 2-track events, 
(97.49 ± 0.03)% for K e3 events). As a check, the analysis was repeated for 
the Ql * 2trk trigger alone, which was measured to be equally efficient for 
all events. The relative uncertainty due to different trigger efficiencies is very 
small: AR/R = 0.05 %. 

In about 5 % of the events the drift chambers record multiple hits in one layer 
which lead to an overflow condition. This could be more likely for electrons 
than for minimal ionizing pions or muons. Comparing the results with or 
without cutting on the overflow condition shows that the effect on R is almost 
negligible: AR/R = 0.05 %. 

Using a data set of monochromatic single pions or electrons from a test run, 
it has been checked that the efficiencies to record and reconstruct pions and 
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electrons are equal within 0.05%. 



In order to be independent of potential asymmetries in the setup, about half of 
the data were recorded with positive polarity and half with negative polarity 
of the spectrometer magnet. We analyzed the data separately for both polar- 
ities, but found an almost negligible difference, resulting in an uncertainty of 
AR/R = 0.07%. 

As a further systematic check the analysis was repeated, broadening a num- 
ber of detector resolutions in Monte Carlo. Energy, momentum and vertex 
positions were convoluted with gaussian distributions, the chosen standard 
deviations being half of the experimental resolution. The number of selected 
events changed only by the order of 10~ 5 , proving that the result does not 
depend on resolution effects. 





relative uncertainty [%] 


experimental normalization (energy spectrum) 


0.67 


normalization error from input ratios 


0.16 


E/p cut 


0.05 


trigger efficiency 


0.05 


DCH overflows 


0.05 


magnet polarity 


0.07 



Table 3 

Summary of systematic uncertainties on the ratio R. 

We summarize the systematic uncertainties on R in Table 3. Using the accep- 
tances given in Table 1, the ratios of branching fractions from Eq. (7) to (10) 
and the above evaluation of the systematic uncertainties, we obtain as average 
two track acceptance a 2 r = 0.2412 ± 0.0004, and a systematic uncertainty in 
the ratio a e /a 2 T of 0.68 %. 



4 Results 



The electron identification inefficiency increases the number of K e3 decays by 
0.49%, while background from misidentified K^ 3 and K 3n decays reduces the 
number by 0.58%, leading to a net correction of -5705 events. This gives: 



B(K L -> neu) 6753478/0.2599 , N 

R = all 2-trLk)) = 12592096/0.2412 = 04978 ± °-° 035 ' < U > 
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As mentioned previously, a second independent analysis resulted in a value of 
R which differs by less than 0.001 from the value in Eq. (11) . 



For the branching ratio of the 3n° decay, the current experimental situation is 
unsatisfactory. We use a weighted mean of the PDG2004 value (21.05±0.28) % 

[I] and the recent measurement of the KTeV collaboration, (19.45 ± 0.18) % 

[II] , and obtain (19.92 ± 0.70) %, where the error is enlarged because of the 
poor agreement of the measurements. Therefore the branching ratio for all 
2-track events is B(2T) = (80.56 ± 0.70) % and 



B(e3) = — ^ = R * B ( 2T ) = °- 4010 ± °- 0028 ± 0.0035, (12) 

T(Kl — > all) 

with the first error being the complete experimental error and the second the 
external error from the normalization, to be combined to 



B(e3) = 0.4010 ± 0.0045. (13) 

This measurement depends on three other measurements of ratios of partial 
K e3 decay widths. This dependence is given by: 



A£?(e3) = ( - O.QQo] * 0.077 - ( - 0.309^1 * 0.075 
V r ( e3 ) / V r ( e3 ) / 

- (^7^r- 0.515^* 0.151. (14) 
VT(e3) J 1 ; 



The decay rate of Kl — * vrez/ is obtained by using the Kl lifetime t{K£) = 
(5.15 ±0.04) x 10- 8 s [1]: 



T(K e3 ) =B(e3)/r(K L ) = (7.79 ±0.11) x lO^s" 1 . (15) 



5 Value of \V US \ 



The CKM matrix element \V US \ can be extracted from the K® 3 decay param- 
eters by ref. [7] 



\V US \ 



\ G 2 F M 5 KO S EW I K o ff*- 



(16) 
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Three quantities in this equation are taken from theory. Sew is the short 
distance enhancement factor, I K o is the phase space integral and f+ n is the 
form factor. 

To determine \V US \ we follow the prescription and use the numerical results 
in ref. [7], where a detailed numerical study of the K e3 decays to 0(p & ) in 
chiral pertubation theory with virtual photons and leptons is presented. The 
integrals given therein correspond to the specific prescription to accept only 
those radiative events which have pion and electron energies within the whole 
K e3 Dalitz plot. From a Monte Carlo simulation we obtain this correction to 
be small 



Number of -fCsr-vi events inside Dalitz plot 

— - c - — = 0.99423. 17 

Number of all K e3 ( 7 ) events 



Using equations (15) and (17), S EW = 1.0232, I K o = 0.10339 ± 0.00063 we 
obtain a value for the product of the CKM matrix element \V US \ and the vector 
form factor f+° n , 

|K*|/+(0) = 0.2146 ±0.0016. (18) 

For the vector form factor, different theoretical calculations have been pub- 
lished recently. Chiral models including the corrections to the order p 6 give 
/ + (0) = 0.981 ±0.010 [7], /+(0) = 0.976 ±0.010 [8] and / + (0) = 0.974 ±0.011 
[9], to be compared with the older value /+(0) = 0.961 ±0.010 [6]. Lattice cal- 
culations in the quenched fermion approximation give /+(0) = 0.961 ± 0.009 
[10], but this value does not include electromagnetic corrections. Taking the 
value from [7], which takes into account chiral corrections to the order p 6 , 
isospin corrections and electromagnetic corrections, we obtain the CKM ele- 
ment to be 

\V US \ = 0.2187 ±0.0028. (19) 



The error on \V US \ is dominated by the theoretical uncertainties, the error on 
f+ (>7r alone contributing ±0.0023. 



6 Conclusions 



We have made a direct measurement of the ratio of K® 3 to all K\ decays with 
two charged tracks, 
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Using the current experimental knowledge of the 3tt° branching ratio, this 
leads to a branching ratio B(e3) = 0.4010 ± 0.0045. This exceeds the PDG 
value by (3.3 ± 1.3)%, or 2.5 standard deviations. It leads to |V^ s |/+(0) = 
0.2146 ± 0.0016, in good agreement with the recent KTeV result [11], but 
larger than the PDG value [1]. Inferring the most recent theoretical evaluation 
of /+(0) = 0.981 ± 0.010 [7], the coupling constant comes out to be \V US \ = 
0.2187 ±0.0028, where the dominant uncertainty is theoretical. This is still 2.4 
sigma lower than required by the 3-generation unitarity of the CKM matrix. 
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